Compatible Finite Elements For Glacier Modeling

We describe the first application of two mixed finite element methods to equations glacier evolution under different simplifying assumptions, along with a framework for implicit solution coupled velocity-thickness equations. The method uses Raviart-Thomas elements velocity and piecewise constants thickness is reframing classic staggered-grid difference case unstructured triangular meshes. second Mardal-Tai-Winther exhibits several desirable properties: order convergence near-exact mass conservation, while resolving both membrane shear stresses in steep topography thin ice.